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Abstract: If the voids that we see today in the distribution of galaxies existed at recombination, 
they will leave an imprint on the cosmic microwave background (CMB). On the other hand, if these 
voids formed much later, their effect on the CMB will be negligible and will not be observed with the 
current generation of experiments. In this paper presented at the 2004 Annual Scientific Meeting of 
the Astronomical Society of Australia, we discuss our ongoing investigations into voids of primordial 
origin. We show that if voids in the cold dark matter distribution existed at the epoch of decoupling, 
they could contribute significantly to the apparent rise in CMB power on small scales detected by the 
Cosmic Background Imager (CBI) Deep Field. Here we present our improved method for predicting 
the effects of primordial voids on the CMB in which we treat a void as an external source in the cold 
dark matter (CDM) distribution employing a Boltzmann solver. Our improved predictions include 
the effects of a cosmological constant (A) and acoustic oscillations generated by voids at early times. 
We find that models with relatively large voids on the last scattering surface predict too much CMB 
power in an Einstein-de Sitter background cosmology but could be consistent with the current CMB 
observations in a ACDM universe. 
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1 Introduction 

Analyses of surveys such as the 2-degree Field Galaxy 
Redshift Survey and the Sloan Digital Sky Survey indi¬ 
cate large volumes of relatively empt y space, or voi ds, 
in the distr ibution of galaxies (e.g. iHovle fc Vogelevl 
2002, I2004T) . In the standard hierarchical model of 
structure formation, gravitational clustering is respon¬ 
sible for empt ying these voids of mass and galaxies 
( Peeble s: Il98flh . However, standard cold dark mat¬ 
ter (CDM) simulations predict significant clumps of 
substructure within voids that that should be capa¬ 
ble of developing into observable bound objects, such 
as dwarf galaxies, that are not apparent in the surveys 

dPekel fc Silkll986HHoffman et alll992tlShandarin et alj 

120044. [Peebles (12001.1 gives an in-depth discussion of 

the contradictions of this prediction with observation. 

He believes that the inability of the CDM models to 
produce the observed voids constitutes a true crisis for 
these models. 

It has been argued that the apparent suppression 
of void object formation may be addressed when more 
sophisticated models for the complex ef fects of baryons 
on ga laxy formation are developed te.g. lvan den Bosh et all 
l200dh . Others believe that altern atives to CDM, such 
as warm dark matter (WDM) (|Hoga n fc Da.lcantonl 

l200 d: ISommer-Larsen fc Dolgovi2fl0lll or self-interacting 

dark matter dSpemer^^Steinh^^20(iflh , may provide 
a solution. The higher temperature of WDM particles 
causes them to move faster than particles of CDM. 


This motion might enable WDM particles to resist 
congregat ing to form obs e rvable galaxies within voids. 
However, iBarkana et all 1 2(111 1 i note that the mate¬ 
rial’s resistance to clumping might delay the epoch of 
quasar formation and the model is ruled out if reion¬ 
isation occurred before z ~ 10 as indicated by the 
recen t WMAP cosmic m icrowave background (CMB) 
data dSoereel et alJl200^4 . 

Self-interacting dark matter may explain the rela¬ 
tive dearth of dwarf galaxies. If there were interpar¬ 
ticle collisions, the halo of dark matter surrounding a 
large galaxy would interact with the halos of nearby 
dwarf galaxies, stripping the dwarves of their gas and 
stars more rapidly than in the standard CDM theory. 
However, Chandra X-ray Observatory data indicates 
that the density of the dark ma tter is greater the close r 
it is to the center of the cluster (lArahadiis et a,ll2002fl . 
in contradiction with the simplest self-interacting dark 
matter models which predict particle collisions to oc¬ 
cur most frequently in more populated regions and act 
to spread out galactic cores and reduce their density. 

The standard CDM model of structure formation 
remains extremely successful at matching the observa¬ 
tions on large scales such as the properties of galaxy 
clusters, the statistics of the Lyman-a forest and the 
temperature anisotropy of the CMB. However, the lack 
of observable dwarf galaxies within voids as predicted 
by the simulations can not be overlooked. Addition¬ 
ally, recent deep field obse rvations from the Cosmic 
Background Imager (CBI) iMason et alJ ll2003f) show 
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excess power in the CMB over the standard model on 
small angular scales, i > 2000. 

It may be possible to explain both these observa¬ 
tions by postulating the presence of a network of voids 
originating from primordial b ubbles of true vacuum 
that nucleated during inflation ltLall99li:fLiddle fc Wandsl 
I1991IL In this scenario, the first bubbles to nucleate 
are stretched by the remaining inflation to cosmolog¬ 
ical scales. The largest voids may have had insuffi¬ 
cient time to thermalise before decoupling and may 
persist to the present day. Such primordial voids are 
predicted to produce a me asurab l e cont rib ution to the 
CMB on small-scales fsee ICriffiths et al. (2003) and 
references therein). If voids in the matter distribu¬ 
tion were present prior to the onset of gravitational 
collapse, the formation of bound objects within voids 
would be suppressed. 

_ In this paper, we discuss the results of lOriffiths et al.l 

( 2003 1 in which we develop a general method to ap¬ 
proximate the void contribution to the CMB that al¬ 
lows the creation of maps and enables us to consider an 
arbitrary distribution of void sizes. We show that if the 
voids that we see in galaxy surveys today existed at the 
epoch of decoupling, they could contribute significant 
additional power to the CMB angular p ower spectrum 
between 2000 < t < 3000. Further to ICriffiths et all 
( 2003 1. we show how we improve our predictions. We 
include the effects of the cosmological constant as well 
as oscillations in the matter-radiation fluid that may 
be generated by primord ial voids on scales up to the 
sound horizon ifBaccigalnni fc Perrottal2000ll and present 
our method for constraining the parameters of this in¬ 
flationary model of void production. 

2 First Approximation 

2.1 Void parameters 

In ICriffiths et all ll2003f) we model the voids seen to¬ 
day as spherical underdensities of 5p/p = — 1. As a 
primordial void bubble nucleates and is stretched by 
inflation it will sweep up the surrounding energy den¬ 
sity to form a thin boundary wall. At the end of the 
inflationary epoch, any voids that persist will be fully 
compensated. This means that the overall cosmology 
is that of the background universe, since a compen¬ 
sated void does not distort space-time outside of itself. 
This is an extremely important property, as it allows 
us to place many voids into a universe without the 
worry that they might influence each other. As a first 
approximation we take the background u niverse to be 
an Ein stei n-de Sitter (EdS) co smology. iMaeda fc Satol 
lll983fl and lBertschingerl Il995h use conservation of mo¬ 
mentum and energy respectively to show that com¬ 
pensated voids in an EdS background cosmology will 
increase in radius r v between the onset of the gravi¬ 
tational collapse of matter at equality and the present 
day such that 

rv{v) oc rf , (1) 

with (3 ~ 0.39 and where r/ is conformal time. 

Motivated by the extended inflationary scenario 
(lLa fc Steinhardtll989l . see Kolb 1991 for a review), we 
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assume a power-law distribution of bubble sizes greater 
than a given radius r of the form, 

Nb{> r) oc r~ a . (2) 

Typically, extended inflation is implemented within 
the fr amework of a Jordan-Brans-Dicke theory llBrans fc Dicke! 
Il96lh . In this case, the exponent a is directly related 
to the gravitational coupling u> of the scalar field that 
drives inflation, 


Values of lo > 3500 are required by solar system ex¬ 
periments (Will 2001). We take the limit of large u>, 
lea ding to a spectrum of vo id sizes with a = 3. 

iLiddle & Wand3 show that a distribution 

of void sizes that predicts the existence of arbitrar¬ 
ily large inflationary voids will cause significant effects 
on the CMB that contradict observations. We assume 
that the mechanism creating the voids imposes an up¬ 
per cut-off on the size distribution. A possible mecha¬ 
nism for this cut-off could be that the tunneling proba¬ 
bility of inflationary bubbles is modulated through the 
coupling to another field. 

As well as avoiding the well known problems asso¬ 
ciated with arbitrarily large voids existing at the epoch 
of decoupling, this assumption allows us to match the 
observed upper limit on void sizes from the galaxy 
redshift surveys. Since the results of void finder al¬ 
gorithms in the literature include those obtained from 
incomplete galaxy surveys, we choose to take the av¬ 
erage maxi mum void size that is found, r max = 25 
Mpc /h ie.g. iHovle fc Vogelevll2003 . 120041 . An analy¬ 
sis is also in progress in which this parameter is kept 
free so that it may be constrained by the observations. 
The minimal present void size is also chosen to agree 
with redshift surveys, r m i n = 10 Mpc//i. 

We normalise the distribution by requiring that the 
total number of voids satisfies the observed fraction of 
the universe filled with voids today, F v . Redshift sur¬ 
veys indicate that approximately 40% of the fractional 
volume of the universe is in t he form of voids of un- 
derde nsity 8p/p < —0.9 (e.g. IHovle fc Voeelevl 120021 . 
I2004I1 . ie. F v ~ 0.4. The positions of the voids are 
then assigned randomly, making sure that they do not 
overlap. In order to speed up this process, we consider 
only a 10° cone. This limits our analysis to f > 100, 
which is satisfactory for our purpose since the main 
contribution from voids is on much smaller scales. 

2.2 Stepping through the void net¬ 
work 

Refer to iGriffiths et all (l2003Tl for a more detailed de¬ 
scription of our methodology. We ray trace photon 
paths from us to the last scattering surface (LSS) for 
the 10° cone in steps of 1’. Each void in the present 
day distribution that is intersected by the photon path 
is evolved back in time according to equation 0 to de¬ 
termine whether the photon encounters the void. 

If a photon intersects a void between us and the 
LSS, we compute the Rees-Sciama (RS) (1968) effect 
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Figure 1: The map of the temperature fluctuations 
on the surface of last scattering from the fiducial 
EdS void model considered. The parameters of 
this model are given by a = 3, r max = 25 Mpc /h 
and F v = 0.4. 


due to the deviation in the redshift of the photon as 
it passes through the expanding void and the lensing 
effect due to the deviation in its path. If a photon 
intersects a void on the LSS, we calculate the Sachs- 
Wolfe (SW) (1967) effect due to the photon originating 
from within the underdensity. We take into account 
the finite thickness of the LSS, which suppresses the 
SW effect for small voids, by averaging the contribu¬ 
tion from a number of photons originating from a LSS 
of mean redshift 1100 and standard deviation in red- 
shift 80. We also calculate the partial RS effect that 
arises due to the expansion of the void on the LSS as 
the photon leaves it. 

Once the photon has reached the last scattering 
surface, we know the variation of its temperature as 
well as its position on the LSS and can create a tem¬ 
perature map (figure 0. We then use a flat sky ap- 

tain the Ci spectrum of the anisotropies (see figure 
0 We point out that primordial void parameters are 
still poorly constrained by both observation and the¬ 
ory. The bottom panel of figure 0 shows a few further 
example models. 

For a power-law size distribution (as motivated by 
the inflationary scenario), large voids become rarer as 
a is increased. Therefore, since void analyses of red¬ 
shift surveys only sample a fraction of the volume of 
the universe, there may exist voids of larger r max than 
currently observed. Models with high r max voids in an 
EdS background cosmology tend to predict too much 
power on scales t ss 1000. However, if we take in- 
flationary model s _with a > 6, as motivated by eg. 
lOccionero fc Amendolal t !994f) . then the peak moves 
to larger £ and the total power drops. The filling frac¬ 
tion mainly adjusts the overall power. 



Figure 2: Top: The CMB anisotropies produced 
by the fiducial EdS void model (solid line) com¬ 
pared to the primary ACDM CMB anisotropies 
(dashed-triple-dotted). Also plotted are the sum 
of primary and void contributions (dotted) as well 
as the fluctuations induced purely by voids on the 
last scattering surface (dashed) and by those be¬ 
tween last scattering and today (dashed-dotted). 
We show the “standard” cosmological concordance 
model: of course a combined analysis of primary 
and void-induced fluctuations would select a dif¬ 
ferent cosmology for the primary contribution. 
Bottom: Example models depicting a range of 
void contributions to the CMB fluctuations. The 
models plotted are a = 3, r max = 25 Mpc /h and 
F v = 0.4 (solid line), a = 3, r max = 40 Mpc /h and 
F v = 0.4 (dotted), a = 3, r max = 40 Mpc /h and 
F v = 0.2 (dashed) and a = 6, r max = 40 Mpc /h 
and F v = 0.4 (dashed-dotted). 


3 Improved Prediction 

3.1 Void Evolution in a ACDM back¬ 
ground cosmology 

The addition of a cosmological constant (A) to the cos¬ 
mology is expected to slow down the conformal evolu¬ 
tion of the voids at late times. We test this hypoth¬ 
esis by simulating the comoving evolution of a single 
void of radius 25 Mpc/ft starting at z = 1000 in a 
(100 Mpc/ft) 3 box containing 64 3 particles. We find 
that the void size evolution in a ACDM background 
deviates from that of the EdS scaling solution at late 
times as expected. However, the final radius is un¬ 
derestimated by less than 2% (see figure 0. The EdS 
void scaling relationship given by equation 0 with 
j3 = 0.39 is therefore a good approximation for a void 
in a ACDM background. 

Since a ACDM universe evolves for substantially 
longer, we would expect the voids that we see today to 
appear smaller on the last scattering surface for this 
cosmology than in our first approximation EdS back¬ 
ground. We would therefore predict that voids in a 
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Figure 3: The simulated size evolution for a 25 
Mpc /h void in a ACDM background cosmology 
(triangles) compared with the scaling relationship 
(equation[T]with P = 0.39) for an EdS background 
(crosses). The two evolutions deviate by less than 
2% at ?y = 0.046, the current age of the universe. 


ACDM universe will have a more suppressed effect on 
the CMB than we have first approximated. This can 
be modelled by taking the horizon size of our cone 
of voids to be that of a ACDM cosmology. Figure 0 
compares the predicted power from a uniform distri¬ 
bution of equally sized voids in EdS and ACDM back¬ 
ground cosmologies. The overall contribution to the 
CMB from the voids is suppressed as expected and the 
peak is moved to smaller angular scales since the voids 
now appear smaller on the last scattering surface. 

3.2 Modelling acoustic oscillations 

So far we have studied voids which do not interact with 
their surroundings. Specifically, they do not perturb 
the matter-radiation fluid themselves. While this is 
a reasonable approximation for the RS effect after re- 


Figure 4: The predicted CMB power from a uni¬ 
form distribution of equally sized voids in EdS 
(dotted) and ACDM (solid) background universes. 
The voids have a size of 40 Mpc /h and a filling 
fraction of 40 %. The void contribution is com¬ 
puted by ray tracing a void network. The effect of 
A is to suppress the effect of the voids and move 
the peak in the power spectrum to smaller angu¬ 
lar scales since they will appear smaller on the last 
scattering surface. 


ternal source (the voids) and the matter-radiation fluid 
is via the metric perturbations <I> and ’3/ and, since they 
are small during the times of interest, they can hence 
be handled within a perturbative framework. 

As opposed to defects, which are entirely external 
sources, the voids represent a perturbation of the dark 
matter. By construction, the dark matter is in a stable 
configuration and does not perturb itself. It is there¬ 
fore influenced by the presence of the void. Also, since 
there is no dark matter inside the void, there can be 
no flow outwards. We have chosen to model these ef¬ 
fects by making the approximation that the dark mat¬ 
ter is completely decoupled from the rest of the uni¬ 


combination, it has been pointed out iBacci^hmifc PerrottJ erse - Therefore, the cold dark matter acts only as an 


|2000|) that this is not true before last scattering. A 
void could therefore set up oscillations in the matter- 
radiation fluid on scales up to the sound horizon, which 
might be clearly visible in the power spectrum. 

The standard way to include the full behaviour 
of matter and radiation perturbations is to solve the 
Boltzmann equation numerically. But it is not possible 
to model a void self-consistently inside the Boltzmann 
solver. The dark matter density contrast inside the 
void becomes very large (approaching —1 at late times) 
and is clearly non-perturbative. A similar problem was 
encountered a few years ago w hile investigating topo¬ 
logical defects fsee lDurrer et all ll2002f) and references 
therein). The solution in this case was to model the 
defects as external sources. This can be done for the 
voids as well. Any back reaction from the perturbation 
on the voids should only be second order so can be ne¬ 
glected. Furthermore, the connection between the ex¬ 


external seed where the perturbations are concerned 
(of course it is taken into account for the background 
quantities). 

We assume that perturbations driving away from 
the void configuration are only due to the back-reaction 
from the induced perturbations in the other fluids. As 
these are second order effects, we can neglect them. 
We could take them into account by coupling the dark 
matter to the fluid perturbations alone. However, as 
expected, the effect on the resulting power spectrum 
is minor, around 1% for the scales of interest which is 
comparable to the accuracy of the Boltzmann codes. 
We prefer to use first order perturbation theory con¬ 
sistently and turn the dark matter perturbations off 
completely. 

We assume that the photons and baryons have 
flown relativistically back into the void, while the cold 
dark matter remained frozen in the wall. We therefore 
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Figure 5: The behaviour of the perturbations in 
real space at last scattering (LS) from a single 
void. The void, which has a radius of 40 Mpc to¬ 
day, has a comoving size of 9 Mpc at LS. The black 
solid curve is the input perturbation $s(x, Vls)/2, 
while the red dotted curve is the total metric 
perturbation <E>/2. The blue small dashed curve 
shows —2'F = 2$. The cyan large dashed curve is 
the temperature perturbation from the photons, 
D g /4. The baryon velocity is shown as the ma¬ 
genta dot-dashed curve. The green small dashed- 
large dashed curve is D g / 4 — 2$, which is very 
similar to $ 5 / 2 . 


choose a density contrast 


S 


_ Pc _ 

P7 + Pu + pb + pc 


Me 

flrel/n + -f- f lb 


(4) 


In the phenomenological model investigated here, the 
voids emerge slowly from the isotropic background of 
relativistic particles during inflation. Therefore, equa¬ 
tion [1] also describes the void density contrast at high 
redshifts. A more violent formation history can lead 
to stronger acoustic waves and cha nge our results. 

Fol lowing the form alism of lKunz fc DurreH f l997tl. 
iDurrer fc Sakellariadoui (Il997f) an d lPnrrer et all ill 
the source is found to be 



Figure 6 : A “zoom” onto the sound waves gen¬ 
erated by the void. Notice that the x axis now 
extends beyond the sound horizon (at about 111 
Mpc) and that the y axis has been magnified by 
roughly a factor 1000 compared to fig. El The lines 
are the same as in that figure. 


of the void. With our form of 5(r]), given by equa¬ 
tion we see that the amplitude of the sound waves 
is only about 1% of the amplitude of the temperature 
perturbations inside the void. For smaller voids, the 
relative amplitude will be even less. This seems quite 
unimportant, indeed, the oscillations are not even vis¬ 
ible in figure El However, the void itself covers a sur¬ 
face of only about 9 2 7r(Mpc /h) 2 on the LSS, the sound 
horizon in our flat matter dominated universe of the 
example case is at lllMpc//i and the sound waves ex¬ 
tend out to about 100Mpc/h, ten times further than 
the size of the void. Therefore, the total power in the 
fluctuations is about 100 times larger than expected 
from their amplitude and is comparable to the the 
power from the void predicted by our first approxi¬ 
mation method. Indeed, we see oscillations appear on 
large angular scales in the angular power spectrum (see 
figure 0. 

3.3 Correlation function approach 


$(77, k) 


An 




5f(k,r v (ri)) (5) 

.( 6 ) 


(3 — ( kr ) 2 ) sin(fcr) — 3kr cos(kr) 


k s 


We inse rt t his sourc e into a modified version of CM- 
BEasy (iDora nl l2003li . We then write out a snapshot 
of the fluids at the time of last scattering and Fourier 
transform them back. 

Figure 0 shows a close up of the perturbations 
around the sound horizon. These are due to oscilla¬ 
tions set up during the formation of the void and are 
therefore strongly dependent on the formation history 


In order to calculate the angular power spectrum of 
the CMB, we can again use a technique developed for 
studying topological defe cts, n amely the use of unequal 
time correlators iTurokl |1996j). The modified Boltz¬ 
mann equation described above is just a system of lin¬ 
ear equations, 

DA'=5. (7) 

Here V is a time dependent linear differential opera¬ 
tor, A' is the vector of perturbation variables and S 
the source term. For one single void, the source is de¬ 
terministic. But in general we are interested in a void 
network with a random distribution of both the po¬ 
sition and sizes of the voids. In this case the source 
becomes a random variable. 
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For any given initial conditions we can formally 
solve the Boltzmann equation by means of a Greens 
function Q, 

/ to 

dtGj(to,t, k)S(t,k). (8) 

The power spectrum is then a quadratic expectation 
value of the form 

C t ~ <X,-(to,k)X,*(to,k)) = 

f*° di0 i (to,t,k) f t 0 dt'gi(t 0 ,t',k)(S(t,k)S*(t',k)) , 

(9) 

and we need to study the unequal time correlation 
function (UTC) of the source in more detail, 

(®(i,k)**(t',k)). (10) 

As the $ of the different voids don’t overlap in real 
space, we can just add them up to create the total <f> 
of the void network, 


5> 

j 

(11) 

5> ( ^(x-cd. 

(12) 


3 


The sum runs over all voids, Rj is today’s radius of 
void j and c j is the centre of void j. Going into Fourier 


space we find 



$(f,k) = 

J dx<I>(f, x)e lkx 

(13) 

= 

'y ' J dx^ Rj \t,x — c_,)e' kx 

j 

(14) 

= 

j 

(15) 

The function k) is just the Fourier space source 


for a single void at the origin with radius Rj at time 
t as given in equation ®. The different spatial po¬ 
sitions of the voids are now encoded only in an addi¬ 
tional phase factor. The UTC required for computing 
the power spectrum is 

Y $ ( Rj) (t , kk)e ik(Cj-Cl) . (16) 

3,1 

Let us first discuss the special case of equal-sized 
voids, Rj = R\/j. In this case we can move the metric 
perturbation $ outside the sum and the UTC can be 
separated into two parts, 

(S(f, k)5*(t', k)> . = 4> (fl) (t, k)/(k) <F (J?) ‘(t', k)/(k) 

(17) 

where the <f>^ are again the single void metric per¬ 
turbations of equation ilt and 


/(k) = Y eiHcj ~ Cl) ( 18 ) 
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The contribution to / can be naturally separated into 
two different parts, namely the contribution with i = 
j which amounts to VAvoid, the square root of the 
number of voids in the universe, and the contribution 
with i ^ j which is a superposition of oscillations with 
different frequencies, given by the distance between the 
voids. These oscillations will in general tend to cancel 
out so that we are left with /( k) ss y/N vo id . 

The fact that the UTC separates into two parts is 
crucial, since it permits us to integrate independently 
over t and t' in equation Each integral amounts 
to the same as solving the Boltzmann equation with 
the source 

S(t,k) = <F (i?) (t,k)/(k) (19) 

and since they are the same, the angular power spec¬ 
trum is just the one computed with this source. 

By setting /(k) ss VAvoid as discussed above, and 
writing C\ R ^ for the power spectrum which we get with 
a single void as a source, we find 

Cl M Avoid(20) 

and not Ce « N^ oid C^ as might have been guessed 
naively, the power spectrum being quadratic in 8T/T. 

In the general case of equation II61 we cannot sepa¬ 
rate the UTC into two parts. But now the oscillations 
in the exp{ik(cj — c;)} term are being helped to av¬ 
erage out by the now incoherent oscillations from the 
<J>(A)) with different radii (as long as the void radii and 
the void positions are not correlated). We therefore ex¬ 
pect the main contribution to the source to come from 
i = j, leaving us with 

^$ ( ^ ) (t,k)$ ( ^ ) (t , ,k) (21) 

3 

as the source. This is a sum of a source which separates 
into two equal parts. Again writing for the single¬ 
void power spectrum, we find 

C e =Y c e Rj) - (22) 

3 

This can be summed up using an actual realisation of 
a void distribution. But since we are mainly interested 
in the average result for the theoretical distribution, we 
can use this distribution directly. Writing p(R)dR. for 
the number of voids with radii between R and R + dR 
the total power spectrum becomes 

C e = J dRp{R)c[ R) . (23) 

3.4 Boltzmann imprint of primordial 
voids on the power spectrum 

Figure |7| compares the void Ce from ray tracing and 
the Boltzmann approach for a uniform distribution of 
40 Mpc /h sized voids. The overall effect from the first 
approximation ray tracing method appears to be sup¬ 
pressed using the improved Boltzmann solution. Evi¬ 
dence of acoustic oscillations generated in the photon- 
baryon fluid can also be seen on large angular scales. 
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Figure 7: We compare the Cg from ray tracing 
and the Boltzmann approach. The dotted curve 
shows the total ST for a ACDM model with voids 
of a size of 40 Mpc /h and a filling fraction of 40 
%, computed by ray tracing a void network. The 
solid curve shows the same but computed with a 
Boltzmann code. The small and large dashed lines 
show the RS effect due to alone for both meth¬ 
ods and with the same approximations, these two 
curves should coincide. 

The RS effect due to <f>„ alone, using the same approx¬ 
imations in both calculations, is in agreement for both 
methods as expected. 

The suppression that is evident using the Boltz¬ 
mann solution over the ray tracing approach is due 
to a number of factors. In the ray tracing method 
we assume that the density contrast 5 = — 1 and we 
take (a'/o) 2 to be 4/r; 2 . Close to radiation domina¬ 
tion (ie. at early times, just after last scattering, when 
we get the biggest effect) both these approximations 
act to increase the ray tracing result with respect to 
the more accurate Boltzmann prediction, together by 
more than a factor of 2. The ray tracing method also 
neglects the Doppler contribution from the baryon ve¬ 
locity, which suppresses the result further. Finally, the 
finite thickness of the LSS and Silk damping were not 
fully taken into account in our first approximation ray 
tracing method. 

The parameter r max is currently fixed using a num¬ 
ber of void analyses of survey s, including the 2-degre e 
Field Galaxy Redsliift survey (iHovle fc Vogelevi2004fl . 
As discussed in subsection l2.2l these surveys only sam¬ 
ple a fraction of the volume of the universe. Therefore, 
the existence of voids of larger r max than currently ob¬ 
served can not be ruled out. Models with a = 3 and 
high r max voids in an EdS background cosmology tend 
to predict too much power on scales £ « 1000. Our 
results show that this is less of a problem for a ACDM 
universe. Furthermore, since the Boltzmann solution 
predicts the CMB power from voids to be lower than 
implied by the ray tracing estimation, relatively large 
voids on the last scattering surface may be consistent 
with the current CMB data, even for void distributions 



Figure 8: Top: The CMB anisotropies produced 
by a ACDM void model (solid line: a = 3, 
r max = 55 Mpc /h and F v = 0.4) compared to 
the primary ACDM anisotropies (dotted). Also 
plotted are the sum of primary and void contri¬ 
butions (dashed). We show the “standard” cos¬ 
mological concordance model. Again a combined 
analysis of primary and void-induced fluctuations 
would select a different cosmology for the primary 
contribution. Bottom: The large dashed curve 
shows the fiducial EdS model of figure 0 (a = 3, 
r max = 25 Mpc /h and F v = 0.4) using our ray 
tracing method. The same model is shown for 
our Boltzmann solution with EdS (small dashed) 
and ACDM (dotted) background universes. The 
solid bold line is the same ACDM void model for 
r max = 55 Mpc//i maximally sized voids as in the 
panel above. 

with low values of a (see figure 0 ■ 

4 Conclusions 

The cosmic microwave background is an excellent tool 
for probing the distribution of matter from last scat¬ 
tering until today. In the case of voids, the strongest 
signal stems from objects at very high redshifts, espe¬ 
cially from those already present at decoupling. We 
discuss in this Paper our improved method to investi¬ 
gate the imprint of a distribution of primordial, spher¬ 
ical and compensated voids, which could for example 
be generated by a phase transition during inflation. 

We show that the signature of a power-law distri¬ 
bution of such voids, that is compatible with redshift 
sureys, contributes additional power on small angular 
scales. The amount of power generated depends on 
the parameters that describe the void size distribution 
as well as the background cosmology. The overall sig¬ 
nal predicted by our first approximation ray tracing 
method is somewhat suppressed using our improved 
Boltzmann solution. Thus, primordial void formation 
models that produce relatively large voids on the last- 
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scattering surface may not be ruled out by the current 
CMB observations. 

Experiments such as the CBI are able to directly 
probe small angular scales and constrain void param¬ 
eters. We will present a Markov Chain Monte Carlo 
constraints analysis of a wide range of void models in 
a future paper. We will further constrain models that 
are comp atible with CMB ob servations using cluster 
evolution llMathis et al.ll2004f) and also investigate the 
non-Gaussian signal of void models that are compati¬ 
ble with the observations. 

Other sources are also expected to contribute at 
high l. Probably the strongest of these is the ther¬ 
mal Sunyaev-Zel’dovich (SZ) effect. Since the ther¬ 
mal SZ effect is strongly frequency-dependent, experi¬ 
ments which work at about 30 GHz (like CBI) will see a 
stro nger signal than those working at higher frequen¬ 
cies llAgha.nim et a,l . ssa). Hence a multi-frequency 
approach should be able to easily disentangle the con¬ 
tribution of voids from the SZ effect. Unfortunately, 
it seems to be difficult at present to predict the pre¬ 
cise level of the SZ contribution, since d ifferent groups 
are re porting different results (see eg. ISnringel et a.lJ 
( 200 if) for a compilation). Future multi-frequency, high 
resolution and high signal-to-noise maps should be able 
to significantly constrain the contribution of primor¬ 
dial voids to the high l CMB power spectrum. Ad¬ 
ditionally, deep galaxy redshift surveys and measure¬ 
ments of the distribution of matter in the Ly-ct forest 
will be able to directly explore the presence of voids in 
the baryonic matter distribution at low redshifts. 
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